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Wave-front tracking for a quasi-linear scalar conservation law with 
hysteresis 


Fabio Bagagiolo*, Stefan Moreti** 


Department of Mathematics, University of Trento, Italy 


Abstract 


In this article we deal with the Cauchy problem for the quasi-linear scalar conservation law 
ut + F(u); + Ue = 0, 


where F is a specific hysteresis operator, namely the Play operator. Hysteresis models a 
rate-independent memory relationship between the input u and its output. Its presence in 
the partial differential equation gives a particular non-local feature to the latter allowing us 
to capture phenomena that may arise in some application fields. Riemann problems and 
the interactions between shock lines are studied and via the so-called Wave-Front Tracking 
method a solution to the Cauchy problem with bounded variation initial data is constructed. 
The solution found satisfies an entropy-like condition, making it the unique solution in the 
class of entropy admissible ones. 


Keywords: Hysteresis, Conservation Laws, Weak Formulation, Entropy Solution, 
Wave-Front Tracking 
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Introduction 


In this work we deal with the Cauchy problem for a conservation law with hysteresis 
nonlinearities as follows: 


Ut + Wi + Uz = O on R x [0, T), 

w(x, t) = Flulz,-), wols) t) Vte(0,T), ae. x €R, (1) 
u(x, 0) = 1g) in R, 

w(x, 0) = wo(z) in R. 
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where: F is a so-called hysteresis operator, namely representing a memory-dependent input- 
output relationship between the pair of scalar functions (t +> u(z,t),t => w(z,t)), one for 
almost every x; uo is the initial datum for the solution u; w? is a suitable space-dependent 
function for the initial values of the output w. The memory dependence represented by F is 
rate-independent, which is the main characterization of the hysteresis phenomena. We are 
going to consider the case where F is the so-called Play operator. 

The presence of the hysteretic term w gives (1) a particular non-local feature which, 
up to the knowledge of the authors, for this kind of equation was not investigated before 
in the framework of characteristics, wave-front tracking and the corresponding limit pro- 
cedure. More precisely, due to the rate-independent memory, the functional input-output 
relationship u œ> w = Flu, w?] is highly non-linear and non-differentiable and this fact leads 
to abrupt changes of characteristics after any possible time t when, in dependence on x, the 
pair (u, w) reaches suitable regions of the phase-space u—w. In the case of the Play operator, 
such abrupt changes can be somehow seen as the case of a discontinuity in the derivative of 
the flux function f for a generic conservation law u;+ f(u), = 0. For the Riemann problem, 
(1) is then seen as a conservation law with a piece-wise linear flux function f, where the 
alternation of the values of f’ depends on the hysteretic relationship u +> w and is different 
for different points x. 

Our main goal is the study of a generic Cauchy problem for (1) with initial data uo, wo 
just bounded variation functions (BV). We then approximate uo and wo by piece-wise 
constant functions, solve the corresponding Riemann problems and, suitably adapting the 
wave-front tracking method, we pass to the limit. In doing that we face several non-standard 
problems: 

i) with respect to the classical wave-front tracking method in our case the number of 
discontinuity waves after an interaction may increase, while the total variation of the solution 
remains constant. 

ii) for the Riemann problem the solution u is discontinuous in time and this is in general 
a problem for the hysteresis relation which is usually defined for continuous inputs u; we 
need to extend the definition of the Play operator to time-piece-wise constant inputs in the 
spirit of the so-called regulated functions (see Brokate-Sprekels [3], Krejci-Laurencot [13] 
and Recupero [18]); 

iii) for the Cauchy problem, being the solution only BV, we need a further relaxations 
of the hysteresis relation writing it as a suitable measure-dependent variational integral 
inequality (suitably adapting the one in Visintin [21]). In particular, such variational in- 
equality will be maintained to the limit, giving the existence of a weak solution for the 
Cauchy problem with hysteresis. 

Our main results are the existence of a solution of (1) with BV data and the uniqueness 
in the class of functions satisfying an entropy condition. 

In the article [21], Visintin studies equation (1) with different hysteresis operator (namely, 
the delayed relay). However, the PDE problem is there dealt by a time discretization method, 
without focusing on characteristic curves. In our article instead we use a different approach 
which is more constructive and gives a concrete idea on how the discontinuities behave and 
how the solution evolves in time. 


Hysteresis is a phenomenon often observed in various natural and engineered systems, 
typically characterized by a lag or delay in the response of a system to changes in the 
input. Well-known examples of hysteresis are in the behaviour of ferromagnetic materials, 
stress-strain relationship in plasto-elastic materials and behaviour of thermostats. For a 
comprehensive account for mathematical models for hysteresis and their use in connection 
with PDEs, we refer to Krasnoselskii-Pokrovskii [12] and Visintin [19]. 

Hyperbolic and scalar conservation laws with hysteresis for some specific applied mo- 
tivations were studied in [16], [17], [10], [11], [4], [15], [5], [8]. In particular, the models 
studied by Peszynska and Showalter in [16] and [15] come from applications in transport 
with adsorption in porous media, where hysteresis is a common feature. Kordulova in [11] 
summarizes the results known for such equations and Kopfova in [10] focuses on entropy 
conditions. Corli and Fan in their works [4] and [5] investigate a conservation law with 
hysteresis relation in the flux given by a parametric family of curves which are followed 
subject to the monotonicity of the input and to the trajectory of an auxiliary ODE. The 
model is motivated by application in traffic flow, where hysteresis is due to a delay in change 
of drivers behavior. Fan in the recent article [8] studies the same model as in [4],[5] and 
considers the Wave-Front tracking method, proving an estimate on the total variations of 
the solutions. However, the limit procedure, in order to have existence of a weak solution 
for more general initial data, is not completely performed especially for what concerns the 
passage to the limit into the hysteresis relationship. 

For general basic theory on scalar conservation laws we refer for example to the book 
by Evans [7]. For more specific results, such as the solution of Riemann problems for scalar 
conservation laws with piecewise linear flux and the classical Wave Front Tracking Method, 
we refer to the books by Bressan [2] and by Holden-Risebro [9] or to the seminal paper by 
Dafermos [6]. 

The article is structured as follows. In Section 1 we introduce the mathematical for- 
malization of hysteresis and we define the Play operator F and its suitable extensions. In 
Section 2 we recall some results for scalar conservation laws and in particular we focus on 
the case when the flux is piece-wise linear. In Section 3 we introduce a suitable formulation 
of (1) and in Section 4 we study the corresponding Riemann problem. In Section 5 we prove 
existence of weak solutions for the case of BV initial data, in particular we perform the limit 
in the wave-front-tracking procedure. In Section 6 we give an entropy condition, showing 
that the solution constructed in the previous section is the only entropy solution. 


1. The Play Operator for hysteresis and its extensions 


Figure 1 represents the so-called input-output Play hysteresis relationship between a 
time-dependent scalar input u and a time-dependent scalar output w. Here the amplitude 
a > 0 is fixed, and we denote by £ the strip {(u,v) : |u — w| < a} which is going to 
represent the feasible states of the system. If at certain time t the pair (u(t), w(t)) satisfies 
|u(t) — w(t)| < a, that is if it belongs in the interior of £, and if the input changes in time, 
then the output will not change, until the pair (u,w) will possibly reach one of the two 
boundary lines of L. If w(t) = u(t) — a, that is the pair (u,w) is on the lower boundary 
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Figure 1: Play operator 


of L and if the input increases, then the output will increase together with u; if instead 
u decreases then w stays constant and the pair (u,w) enters the interior of £L. If (u, w) 
belongs to the upper boundary of £, then the behavior is symmetric, reversing the role of 
the monotonicity of u. 

Given an initial state (u(0),w(0)) € £ and the evolution of u we can then trace the 

evolution w. We can notice that the value of w(t) is not determined pointwisely by u(t), 
indeed for fixed u(t) we have more than one possible value w, depending on the past evolution 
of u. Hence w(t), besides its initial value w(0), is determined by the whole history uljoq, i.e. 
there is a memory effect involved. 
In particular this memory effect is rate-independent which means that the relation between 
u and w does not depend on the speed of u(t), i.e. it does not depend on its derivative. 
Notice that this requirements is essential if we want to draw hysteresis relations as in Figure 
1 and it is a general feature of hysteresis phenomena. 

When the input u is in W'1(0,T), for some T > 0, the Play hysteresis relationship 
described above can be characterized by the following variational inequality for almost every 
t € [0,7]: 


ju—wl<a, w(u-w-—v)>0 W,|v| <a. (1.1) 


Conditions (1.1) can be interpreted in the following way. If |u — w| < a then the term 
u—w-—v can have either positive or negative sign depending on v. Since the inequality has 
to hold true for every v then necessarily w’ = 0. Whereas, if for example u = w + a, then 
u—w—v is always non negative for all v, then also w’ has to be non negative. The opposite 
happens when u = w—a. In [19] (see also [12]) it is shown how (1.1) well defines an operator 
on the space of Sobolev functions which, by a density argument, can be uniquely extended to 
an operator which maps continuous inputs to continuous outputs. Hence the Play operator 
is proved to be as an operator continuously (with respect to the uniform convergence) acting 
as 


where 


ji { (u, wo) € ©°([0,T]) x R| (u(0), wo) € ch 


Moreover such an operator satisfies 
i) Causality: V (u1, Wo), (u2, Wo), such that u = us in [0, t] then 
F(u, Wo)|(t) = [F (u2, wo) (2). 


ii) Rate-independence: Y (u,wo) € £L, Y t € [0,T] if s : [0,7] > [0,7] is an increasing 
homeomorphism, then 


F(u o s, wo)|(t) = [F (u, wo)](s(é)). (1.2) 


iii) Semigroup property: Y (u, wo) € L, V ty < tz € [0,T] setting w(ty) := [F (u, wo)](t1) 
then we have that 


[F (u, Wo) ](t2) = [F (ufti + +), w(t) (t2 — t1). (1.3) 


In the following, we will need to consider F as applied to u, with u solution to a Riemann 
problem for a conservation law. Hence we have to extend the definition of F to piece-wise 
constant inputs with a finite number of jumps. 


Definition 1.1. A function g : [0,7] > R is piece-wise constant if we can find 0 = to < 
tı <---<ty =T and g; E€ R such that 


N 
g(t) Y gilet) (1.4) 
i=1 


We denote by PC([0, T]) the space the piece-wise constant functions on [0, T]. 


Let us consider the simple case of a function u : [0,7] —> R consisting of two constant states, 
u_,U 4, separated by discontinuity at t = t*, and let us suppose that w_ < u4. The idea is 
to suitably approximate u with continuous functions u, and then apply the operator F to 
these functions and let £ go to 0. So, for € > 0, let us define 


u O0<t<t*-e, 
uelt) = < f(t) t*—e<t<t +e, (1.5) 
UL t*+e<t<T, 


with 


fe : [t — e, t* +£] — [u-, u+] continuous and increasing and 
fe- e) =u- fe +e) = u4. 
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Figure 2: An explicit example of the operator F applied to piecewise constant u. Here a = 1, u- = wo = 0 
and u+ = 2 consequently w* = 1. In red we highlight the path followed by the couple (us, w) which, after 
the limit procedure, collapses to (u_, wo) for t < t* and (u4, w*) for t > t*. 


We can notice that, by the rate-independence (1.2), for every € > 0 the value of output 
w,(t* +£) does not depend on the choice of the function f; satisfying (1.6) (see also Remark 
1.2). Moreover, the value w.(t* +€) =: w* is constant for every £. Indeed changing € in 
(1.5) is just a time rescaling and hence the output of the play operator does not change in 
t* +e. Hence, trivially, limo we(t* + £) = w*. We can then define the output w for the 
discontinuous input u as 


(1.7) 


w? 0<t<t*, 
w* ee ae eel At 


which of course coincides for t 4 t* with 
w(t) := lim we(t) = lim|F (we, wo)](t). 


Note that the output w is also a piece-wise constant function with at most one jump in t*, 
where both u and w do not need to be defined. 

Then with the above construction it is straightforward to extend F to the piece-wise 
constant function so 


F: PC((0,T]) x R > PC((0,T)). 


Remark 1.2. The monotonicity property required to f as in (1.6) is not indeed necessary, 
due to some suitable properties of the Play hysteresis operator (see [19]). 

Also for piece-wise constant inputs and outputs the semigroup property holds as: V (u, wo), 
V ty < t2 € [0,7] such that w(t1) := [F(u, wo)](t1) is defined, then (1.3) holds as equality 
for piece-wise constant functions. 

The way to extend the Play operator to piecewise continuous inputs, as we sketched 
above, can be seen as a special case of a more general issue concerning the extension to the 
so-called regulated functions (see[3] and [18}). 
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In the sequel of the paper we will also perform a limit of solutions of approximating 
Riemann problems, and hence we will need an extension of the Play operator to even less 
regular inputs, as BV. In view of this fact, here we give a further characterization of F for 
piecewise constant inputs which will inspire a weak formulation for BV inputs (see Remark 
3.2). 


Proposition 1.3. Let wo € R and u,w : [0,T] —> R be piecewise constant functions with 
a finite number of discontinuities. Moreover, denote by ù and w the right-continuous rep- 
resentatives of u and w respectively and suppose u(0) = u(0) and w(0) = w(0). Then the 
following are equivalent: 


1. w = F[u, wo] for almost every t; 
2. |\w —u| <a for every t and 


fa — w)d(Dw) > a|Dw|((0,t)) for every t € (0, T], (1.8) 
(0,t) 


where Dw denotes the measure associated to distributional derivative of w, and |Dw| 
its total variation. 


PROOF. First notice that, since wù is of the following form w = ey wilt t) then its 
distributional derivative (and the distributional derivative of w) can be represented by a 
finite sum of Dirac’s delta. In particular 


N-1 


Dw = So (wins — wi)d4,; (1.9) 


i=1 


and 
Dw|((0.2) = S> poi- (1.10) 
{i | ti<t} 
(1 = 2) By definition of F of course |Ù — ŭ| < a everywhere. Let us now consider t € (0, t4], 
then |Dw|((0,t)) = Dw((0,t)) = 0 hence (1.8) holds trivially. Now if t € (t1,t2] we have 
that 


J ©- 0) ade) = Gt) -= WH) (we = w) 

(0,t) 
and 

a|Dw|((0,t)) = alwe — wl. 
Now since w — w # 0, i.e. w did not remain constant, then the arriving point after the 
jump, (u(t ),w2), must be on the boundary of £ (see Figure 2), hence u(t,) — w(t) = +a 
(recall that by definition lim, p+ u(t) = u(t,) and lim, p+ w(t) = w(t,)). In particular we 
can also write (t1) — w(t.) = sign(w2 — w,)a and so (1.8) holds as an equality. Now we 
can can extend recursively the proof on the whole interval (0, 7). 
ve 


(2 = 1) Consider first t € (0,tı] then we have w = wy, in (0,t) and so the constraint 
ù —u| < a tells us that |u — wı| < a for almost every t € (0,t,) so even if u may change 
values it can not exceed the threshold wı + a which would imply a change in w. Hence the 
hysteresis relation F[u, wo] holds. If instead we consider t € (tı, t2] then (1.8) tells us that 


(u(t1) — W(t1))(w2 — wi) > alwe — w)| 
which, because of |w — ù| < a, is equivalent to 
(lti) — W(t1)) = sign(we — wi)a. 


So we conclude that if a jump in w occurs then we reach the upper or lower boundary of 
L, depending whether the jump of w is increasing or decreasing, and of course, due also 
to the constraint (u,w) € £, u must jump with the same monotonicity sign. Hence an 
approximation us as in (1.5) is easily constructed and the corresponding output ws almost 
everywhere converge to w, showing that w = F(u) (see again Figure 2). 


Remark 1.4. In the proof we saw how (1.8) actually holds as an equality. Sometimes, in 
the sequel, it will be convenient to write it as an inequality, nevertheless, because of the 
necessary condition |w — u| < a the two formulation are equivalent. 


Remark 1.5. When we introduce the operator F in a PDE, we have to define it on functions 
u both depending on a space variable x € R”, for some n, and on time t. In this case, for 
every fixed x € R”, we see u(z,-) as a function of time only, and then we define the output 
as 


w(x, t) := (F(u(2,-)wo(x))](t), aex, Vt, (1.11) 


where the initial state of the output is now a given function depending on zx. 


2. Solution to scalar Riemann problems with piecewise linear flux 


In this section we will briefly recall some results for the Riemann problem for a conser- 
vation law with piecewise linear flux function which will be used in Section 4. Let us deal 
with the Cauchy problem 


(2.1) 


u(x, 0) = u(x) xeER, 


iL +f(u,=0  (z,t) €R x (0, +00), 


with f : R > R continuous and piecewise linear and 


<0 
uols) = n - i 


Ur =N, 


With uo as above, (2.1) is usually called Riemann problem. 
A weak solution for (2.1), even for more general f and ug (at least integrable), is here 
defined. 


Definition 2.1. We say that u : R x [0,+00) — R, is a weak solution of (2.1) if u is 


continuous as a function from [0, +00) into Lj,,(R) and 


+00 +00 +00 


f up: + flu)ġa] dx dt + J reeeo dz = 0, (2.2) 


0 —oo —oo 


for every Ct function ¢: R x [0,-+00) + R with compact support. 


It is known that a weak solution with a jump discontinuity along a curve C parametrized 
by t + s(t) with s derivable, should satisfy the so called Rankine-Hugoniot condition (see 
e.g. [7]): 

f(u_(t)) — f(us(t)) = s'(t)(u- (t) — u(t) (2.3) 
on C, where u_(t),u,(t) are respectively the limits from the left and right parts of the 
domain separated by discontinuity curve C. 

In general for conservation laws there is not uniqueness of the weak solution. The intro- 
duction of the so-called entropy condition allows us to select a unique solution. In particular 
one can require that u satisfies 


[ftw tle + sint- DUA) deat > 0, A 


(oe) 


for every k € R and for every non negative smooth ¢ with compact support. It is proven, 
[14] (see also [2], [9]), that for a scalar conservation law there exists one and only one weak 
solution satisfying (2.4). Such unique solution is usually called the entropy solution. 
For piece-wise constant weak solutions, using also the Rankine-Hugoniot condition (2.3), 
the entropy condition (2.4) simply reduces to (see [2]) 
Flu) = Flus) o Flas) = Flu) o Fu) = Fu) Be 


Ur Ut Ut — UL E U— UL 


for every value u between u_ and u+. 

For the case of piecewise linear flux, the entropy solution of the Riemann problem can be 
explicitly constructed, see [6], [2] and [9]. Such solution is piecewise constant with a finite 
number of values. For the purposes of next sections we sketch here such construction. 

We recall that f : R — R is piecewise linear if f is continuous and there is a finite family 
of disjoint intervals, covering R, where f is affine. Moreover u : R x [0, +co|- R is piecewise 
constant if the domain can be partioned in a finite family of subsets where u is constant. 

In (2.1), let us suppose u; < ur, and consider fe the greatest convex minorant of f on 
the interval [u;, u,], which we recall to be 


fe := sup{h convex in [u ur] | h < F}. (2.6) 


It is easy to see that f, is a piecewise linear function hence the graph is a polygonal curve. Let 
us denote the vertices in the following way: (w, f(uz)),(u1, f(u1)), ---, (Uk, f(ur)), (ur, f(ur)) 
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U2 


UZ 


Ul Ur 


Figure 3: An example of solution with piecewise linear flux 


where u < uy < --- < ur and (u, f(u;)), by construction, are some of the vertices of f. 
Because of the convexity we have 


Flu) — flu) — flu) — Fur) —  - Plur) — Flur) 


| (2.7) 
ui — UI U2 — Ui Ur — Uk 
hence we define the following function, see Figure 3, 
Ul itn SNE < Pay fe) 
u Pua) flu) Rex Flua)= flu) 
ui— u u2—U1 
u(x,t) = (2.8) 
m Ny ei) < E < too. 
Ur—uk 


It can be proved that u € C?([0, +00); Zj,,(IR)) and that, by construction, the Rankine- 
Hugoniot condition holds, hence we have a weak solution (as it is known, it is easy to 
see that a piecewise constant function satisfying the Rankine-Hugoniot condition is a weak 
solution). Moreover for every u € [w;, ui+ı] the entropy condition (2.5) holds as an equality, 
due to the piecewise linearity. When u, > u, we consider instead the least concave majorant, 
which again is a piecewise linear function and by considering the vertices of its graph we 


construct the solution similarly as above. 


3. Weak formulation for the Cauchy problem with Play hysteresis 


In this section we consider the following problem 


Ut + We + Uz = O in Rr, 
w = |F (u, wo) |, 
u(x,0) = uo(z) in 


(3.1) 


A 


w(x,0) = wo(xz) inR, 
10 


where Ry := R x [0, T), and the state-dependent hysteresis operator is defined as in Remark 
1.5, starting from the Play operator (with amplitude a > 0) applied to t +> u(x, t) for almost 
every x fixed. 


Definition 3.1. A couple of functions (u, w) with u, w : Rr —> R, continuous as a functions 
from [0, T) into LŁ.(R), is a weak solution to (3.1) if: 


i) it satisfies the following weak formulation of the PDE 


+00 +00 +00 


f fi ar eae E [uote \ + wo(a))o(a,0) de ZO, (3.2) 


—co 


for every C function ¢ with compact support in Ry; 
ii) for almost every (x,t) € Rr it holds 


lu(x, t) — w(x, t)| < a; (3.3) 


iii) for almost every z, the distributional derivative 2 5, is a measure on Rp (denoted in the 
same way) that satisfies 


zee t)? — uo(x)”) dx + ; [ve — wo(z)°) dz < —a E 


for almost every t € (0, T). 


pAg 


R x (0,t)), (3.4) 


Remark 3.2. Similarly to [21], equation (3.4) can be interpreted as an equivalent formula- 
tion of (1.8), in the case of H+ functions which also depend on the state x. Indeed suppose 
that both u, w are H!(Rr) solution of the PDE (3.1) then (3.4) reads as follows 


+00 +00 
J. I wu + ww dtdx + af [ |w,| dtdx < 0, 
so, by the PDE, 


too rt +œ pt 
A ik —(wu + uzu) + ww dtdx + af f |w,| dtdx < 0. 
—oo 0 —oo 0 


Since for almost every t, u is in H'(R) so 


+00 
f u,u dx = 0 for almost every t, 


OO 


we can conclude that 


+oo t +00 t 
i l (w — uj)w, dtdx + af 1 |w,| dtdx < 0, 
—oo 0 —oo 0 


which is indeed (1.8) extended to space-dependent Ht functions. Finally, as in Section 1 for 
time-dependent piecewise constant functions, and arguing as in [20] and [21], (3.4) can be 
seen to be a suitable weak extension of the Play operator to time-space dependent integrable 
functions. 
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If the weak solution has a jump discontinuity on a curve (s(t), t), e.g. say u_(t) 4 u(t) or 
w_(t) Æ w,(t) then we get the following extended Rankine-Hugoniot condition 


u_(t) — u+ (t) ag 
u_(t) — u,(t) + w_(t)— w(t) (t). (3.5) 


The entropy condition (2.4) is also extended to the following one 


f E E ah E SO, (3.6) 
Rr 


for every @ non negative, smooth and with compact support in R x (0,T) and for every 
couple (k,k) € £. 


4. The Riemann problem with Play hysteresis 


We study the Riemann problem for (3.1) where uo(x), wo(x) consist of two constant states 
separated by a discontinuity in the origin. We recall that, for regular input and output, if 
(u, w) belongs to the interior of £ then ws = 0, whereas when (u, w) is on the boundary of 
L then w; = u, whenever w, # 0. Note that on the boundary, we may have w, = 0 and 
uz # 0 only for negligible times t (because in such a case, the pair would immediately enter 
the interior of £L or move along the boundary). 

By the cases discussed here above, we can rewrite (3.1) as follows 


Ut + Uz = 0 ju—w| <a, 
up + tugs = 0 ju — w| =a, 
w= [F(u,wo)) (4.1) 


| 
u(x, 0) = u(x), 


w(x, 0) = w(x 


Of course the initial data must satisfy |uo(x) — wo(x)| < a for almost all z. Our goal is 
to rewrite (4.1) as a unique conservation law with piece-wise linear flux. We consider two 
different cases, a) |uo(a) — wo(x)| < a and b) |uo(x) — wo(x)| = a. Of course a further case 
c) |ug(x) — wo(x)| < a for some values x and |ug(x) — wo(x)| = a for some others, might by 
lacking at a first glance. However, since we are considering a Riemann problem with waves 
moving towards right (see the positives coefficients for u, in the first two lines of (4.1)) the 
cases a) and 6) are the only significant ones. We refer to Remark 4.3 for further explanation. 


a) |uo(x) — wo(x)| < a ae. x 
Let us consider initial data 


U(x) := f TAR (4.2) 


Figure 4: In the first picture we see that if wp — a < u, then we remain inside the hysteresis region when 
jumping from u, to w; in the second picture, when w, — a > uj, instead we hit the boundary, so we have 
w = 0 initially and then uw; = uz. 


with u; < u, and 


w x <0, 
wolz) = f H (4.3) 


hence in this case it is |u — w| < a and |u, — w,| <a. At the intial time t = 0 the PDE is 
then 


Ut + Ur =0, VrER. (4.4) 


If there was no hysteresis the PDE above is solved by the travelling wave solution u(x,t) = 
uolx — t). Hence, if we fix x € R we are expecting our solution u either to remain equal to 
u, for x < 0 or to decrease from u, to u; after a time t = x, when x > 0. 

We have two subcases. 

al) If a is such that w, — a < u then the jump from u, to u is such that the couple 
(u,w) will always remain in the internal region of hysteresis, and so w, = 0 for all z and 
t. Hence we do not experience the effect of hysteresis and so we get the trivial solution, 
u(x,t) = u(x — t) and w(z,t) = wo(z). 

a2) Suppose instead w,—a > uy. For x < 0, u remains constantly equal to u; and so does 
w. For x > 0, we expect our solution to decrease, so if we look at the hysteresis diagram we 
notice that for w, —a < u < u, we have w, = 0, whereas for u; < u < w, — a we are on the 
boundary and so w; = u;. See Figure 4. 

Let us define the following function 


1 0 r ro j 
Nene =l x <0 or (x > 0 and u, > u >w, -— a) (4.5) 


z x > 0 and u <u < w, -—a. 


As already pointed out, for x < 0, the solution is constantly equal to u, (and the output 
to w), hence we can focus on h defined only when x > 0, which is then only dependent on 
u € [u;, ur]. Moreover, the switching rule (4.5) for h actually encodes the hysteresis behavior 
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u: Wr — W: uy +a 


Figure 5: The solution u and corresponding w 


between u and w. We can then rewrite our initial problem as follows 


uz + A(u)u, = 0, 
0) = uol). ie) 


We can see h as the weak derivative of some piecewise linear flux g: [ui,u,]| + R of the 
following form 


glu) == (4.7) 


1 
yu u SUS Up — 4, 
ute Wp — ASU < ur, 


where the constant c makes the flux Lipschitz continuous. Also note that g is convex. Our 
PDE then reads as the following conservation law with piecewise linear flux (compare with 


(2.1)) 
ut + glu)z = 0, (4.8) 


for which we consider the Riemann problem with datum (4.2). As in Section 2, we find the 
solution to (4.8) 


Ul T < z, 
u(x,t) = 4 w,—a 5<2<1, (4.9) 
Up Te. 


We then compute w(2,t) = |F(u(z,-), wo(x))|(t) that, as explained in Section 1 is, see also 
Figure 5. 


WI] x <0, 
w(x, t)= {u t+a 0<2<i, (4.10) 
Wr i< 


Note that the pair (u, w) as above constructed, satisfies the Rankine-Hugoniot condition 
(3.5). Indeed we have two discontinuities on u: the one from u, to w, — a is such that w 
remains constant so 

Uur — (wr — a 
s'(t)= A a ) =f; 
ur — (w, — a) 
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the one from w, — a to u is such that also w jumps from w, to u; +a hence 


P Wr —a-— u 1 
s(t) = =, 
Ww, —a— u +w — (luta) 2 


We also have a discontinuity on w with slope 0 between w, and u; + a but it still satisfies 
the Rankine-Hugoniot condition; indeed u is continuous implying that s'(t) = 0. Being the 
Rankine-Hugoniot condition satisfied by the piece-wise constant pair (u, w), then the latter 
satisfies the weak formulation of (3.1), that is (3.2). 


Remark 4.1. The argumentation above on the Rankine-Hugoniot conditions, works well 
not only because we have posed w = F(u) but mainly because u solves (4.8) where, as 
already pointed out, hysteresis is somehow encoded in the flux g. 

As said, (u, w) can be seen as a weak solution of the Riemann problem (3.1), (4.2), (4.3). 
Also note that, as we already know, the relation w(z,-) = F(u(z,-),wo(x)), for piece-wise 
constant functions t +> (u(x,t), w(x, t)), is equivalent to (1.8). In the proof of Theorem 5.4, 
we will show that (1.8) for piece-wise constant functions, together with the hypothesis of 
integrability on R x [0, T) (which is not satisfied by the solution of the Riemann problem), 
implies condition (3.4), giving a weak solution in the sense of Definition 3.1. 


Remark 4.2. If the initial datum satisfies u; > up, then we end up with a similar situation 
as above but this time the piecewise linear flux g will be concave, coherently with what said 
in Section 2. The construction of the solution then goes similarly as above. 


b) |uo(x) — wo(x)| =a a.e. x 
Consider the initial states (4.2),(4.3) with u; < u,. We have two subcases. 

b1) Suppose w, = u, + a, that is (ur, wp) belongs to the upper boundary of £. Here we 
expect the solution u to decrease from u, to u for x > 0. Then we can notice that the pair 
(u, w) will always remain on the upper boundary w = u +a, see Figure 7. The solved PDE 
is then 


1 
Ut + ple = 0. (4.11) 


The solution u(x,t) = uo(x — 1/2t), and its output w, see Figure 6 are then 


ee w 0< z, 

u t 29 

u(x,t) = E rea owsa Ee (4.12) 
g ee Ur +a i< = 


b2) Suppose (u,,w,) to be on the lower boundary of £, that is w, = u, — a. In this 
case for x > 0 we have u decreasing from ur to u so the couple (u(z,-), w(x, :)) will move 
internally possibly reaching the upper boundary and following it until u = u, see Figure 9. 
We then have two further subcases. 
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Journal Pre-proof 


Ul w 


Ur ur +a 


Figure 6: The solution u and corresponding w for the subcase b1) 


Figure 7: If we start on the upper boundary and we jump from up to u, then we will always remain on that 
boundary. 
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Ur Ur — a 


Figure 8: The solution u and corresponding w in subcase b2ii) 


b2i) Suppose u > w, — a, that is the left-case in Figure 9. By (4.1), when u = u, or 
u < u < ur the PDE is, respectively, 


ut + —Uzp =O or u+ uzr =Q. (4.13) 


So if we define 


1 
5 U = Ur, 
hlu) := 4 ? (4.14) 
1 Ur <U < ur, 
and as in the case a2) consider a conservation law with g piecewise linear flux, then the 
solution is uo(x — t), that is 
u 
u(x,t) = i i 
Ur 


which generates the constant output w(x, t) = wo(z). 
b2ii) Suppose now u; < w, — a so for x > 0 when u decreases from u, to u; we hit also 
the upper boundary at u = w, — a, see the right-case in Figure 9. Again by (4.1) we have 


<1, 


4.15 
eT (4.15) 


+18 +18 


u = Ur, 


h(u) = Wr —a <u < ur, (4.16) 


NIP = nie 


uU < uU < Wp — a, 


which is the derivative of the flux. Then, by solving as in case a2) and subcase b2i) a 
conservation law with piece-wise linear flux and after computing w = F(w,wo), see Figure 
8, we find the pair 


u TID w x <0, 
ulz, t) = $ w,—a Pete. wļ(z,t)=<4u+a 0<2<ł (4.17) 
, r Qo , , l t 2? i 
Ur >l, Ur — a t>% 


As for the case a), Rankine-Hugoniot conditions can be also derived for case b) and 
moreover similar remarks as Remark 4.1 and Remark 4.2 hold. 
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Figure 9: The two possible cases whether wy — a < uj or wr — a > u. 


Remark 4.3. Since we are dealing with a Riemann problem, the possible case c): |wo(x) — 
wo(x)| < a for some values x and |uo(x) — wo(x)| = a for some others, just reduces to either 
lu — w| < a and |u, —w,| = a or ju — w| = a and |u, —w,| < a. Since for x < 0, ulz, t) = u 
and so w(x,t) = w; (the waves are moving towards right) the only significant study is for 
x > 0. Hence the conditions for the cases a) and b) reduce respectively to ju, — w,| < a and 
|u, — wr| = a, and hence, their study also covers case c). 


Remark 4.4. For the subcases al) and b2i), as said, the solution is just a rigid movement 
of the initial data. Instead for the subcases a2) and b2ii) we easily see that u; < w,—a < u, 
and also u; + a is between w, and wy, since w; € [u; — a,u; + a]. For the subcase b1) again 
holds that u; + a is between w, and w;. These facts are important as they imply that for 
every case Var(u(-,t)) = ur — uy = Var(uo(-)) and Var(w(-,t)) = w, — w, = Var(wo(-)) 
for every t, where Var denotes the total variation on R with respect to the spacial variable. 
The fact that the total variation does not increase in time is important and will give a 
compactness tool to prove the existence to the general Cauchy problem as in Section 5. 


5. The general initial data Cauchy problem 


The goal of this section is to construct weak solution in the sense of Definition 3.1. by 
using the Wave Front Tracking Method. 

As already said, such method consist of approximating initial data with piece-wise con- 
stant function, so we consider first 


N 
i=l 
and 
N 
wo(x) = 5 will 1.20) (5.2) 
i=l 


such that |uo(x) — wo(x)| < a where —00 = £o < £1 < +--+ < Ty = +00. 
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eal £2 £3 £4 


Figure 10: A possible solution generated by the above procedure, where we have interactions at time 7, and 
T2. In black the discontinuity lines of slope 1 of u, in red ones of w with slope 0, and both colors for the 
common lines with slope 1/2. 


In order to solve the Cauchy problem (3.1), (5.1), (5.2), we give the following heuristics 
idea, which differs from a classical one by the presence of the hysteretic term w, see Figure 
10: 


1. We solve N — 1 Riemann problems at time t = 0 centered in the points 71,...,%y_1, 
which generate discontinuity lines for u and w; by the analysis in Section 4, we know 
that the discontinuity lines for u are lines with slopes 1 or 1/2 (in the plane t — x), 
whereas the ones for w with slope 0 or 1/2; moreover the discontinuity lines of slope 
1/2 are always discontinuity lines for both u and w. 

2. Let us denote by 7, the first time that either a discontinuity line of u impinges a 
discontinuity line of w with slope 0 or two discontinuity lines of u impinge themselves. 
The pair (u, w) found in the previous point is a weak solution of (3.1) for t < 7 (i. e. 
(3.2) and w = F(u)); note that 7, > 0 exists since the number of constancy intervals 
is finite; 

3. we consider u(-,7,) and w(-,7,) as new initial conditions, since they are still piecewise 
constant, we can solve again a finite number of Riemann problems at {t = tı}. So we 
extended u and w for small times after mı; 

4. we proceed in way extending at each step the solution for larger times. 


Next lemma formalizes the heuristics above. 
Lemma 5.1. For the Cauchy problem (3.1), (5.1), (5.2), with T = +00, arguing as above, 
we can construct two piece-wise constant functions u,w : R x [0, +20) such that: 


i) The couple (u,w) satisfies the weak formulation of the PDE (3.2) and the relation 
w = F(u, wo) for almost every x. 
ii) For fixed t € [0, +20) we have that 


Var(u(-,t)) < Var(ug(-)), Var(w(-,t)) < Var(wo(-)), (5.3) 
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where Var denotes the total variation on R of the functions with respect to the space 
variable. 
iii) For t,t’ € [0, +20) the following estimates hold 


J > lu(x, t) — ulz, t')| dx < Var(uo(-))|t — t'l (5.4) 
and ee i 
J |w(ax,t) — w(x, t") | dz < zV ar(woC))lt —t'|. (5.5) 


PROOF. Att = 0, we solve N—1 Riemann problems as described previously and we compute 
the couple (u, w) for t smaller then 71. By the results of Section 4, we know that (3.2) holds 
and w = F(u, wo) for every x possibly except 71,...,£y—1. By Remark 4.4, (5.3) holds for 
t € [0,7) as an equality. Inequality (5.4) also holds for t,t’ € [0,7,). Again, a possible proof 
of it differs from a standard case ([6], [2], [9]) by the presence of the hysteresis in the equation. 
We give here an idea of that for the single Riemann problem centered at 0 with values u 
and up as in Section 4. Considering the set R := {(x,s) | 1/2 < z/s <1, 0< s < Tı}, we 
can notice that u(-,t) and u(-,t’) are equal to w at the left of R and to u, at the right of 
R. Internally they may be equal either to w, u, or wı + a (see for example (4.9) or (4.12)), 
here we will denote this internal value with u*. Note that, in any case, u* is in between u 
and u,. It is easy to see that if we only have one discontinuity u_ 4 u,, with slope 1 then 


+00 
f E o T EE N (5.6) 


OO 


If instead we have only one discontinuity with slope 1/2 then 


T t r 
/ lu(x, t) — u(x, t')| dr = ; aS 


5 Ju_ — uy] < fé—e'] Ju — uy. (5.7) 


Since we may have at most 2 discontinuities then (5.4) follows from (5.6) and (5.7), indeed 


+00 
f ju(x,t)— u(z,t’)| dx < |t—t'|(\u, — u*| + Ju* — ul) = |t — t| Jur — ul, 
where |u; — u,| = Var (uo). 

Inequality (5.5) is treated in the same way with the difference that discontinuities of 
w have either slope 0 or 1/2 so the maximum is 1/2. Notice that (5.4) and (5.5) imply 
u, w €E C°((0, Ti); Lisel R)). 

The lemma is proven for t € [0, 71) now we need to define u and w for t = 7 and 
use them as new initial data. As a standard non-hysteresis case, the piece-wise constant 
functions u(-,7,), w(-,71) can be constructed as limits as t —> 7, without increasing the 
total variations (see again [6]). Moreover, for x ¢ {%1,...,ev_1} U {y|(y, Tı) belongs to 
discontinuity lines started at t = 0}, it holds w(x, n) = [F (u(x,-), wo(x)|(m) as u and w 
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are constant in (x,t) for t < 7, sufficiently close to tı, see Figure 10. Also note that the 
excluded values of x are a finite quantity. 

Taking the piece-wise constant functions u(-,7,), w(-,71) as initial conditions at 7, and 
repeating the reasoning above, we are able to find a new instant 72 > 7, and to extend u, w 
for t € [™,72), satisfying (5.3), (5.4) and (5.5) in [0,72). Also note that, by the semigroup 
property (see Remark 1.2), w = F(u, wo) in [0, 72) for every x different from the one excluded 
previously. 

Then we apply such extension procedure recursively, getting a sequence of instants Tı < 
T2 < 73 < .... Next two lemmas imply that the number of such instants (and hence of 
recursive cases) is finite and hence we actually may perform the extension up to whole 
time-line [0, +00). Finally note that at each step we exclude a finite number of x such that 
w # F(u, wo), hence w = F(u, wo) in [0, +00) for almost every zx. 


Lemma 5.2. The number of discontinuity lines of u intersecting each other is finite. 


PROOF. Let us denote by w,...,uy and wı,...,wg the initial values of u and w re- 
spectively. At every step-instant 7, by solving the Riemann problems, we end up with 
u taking values in the sets {initial values of u} and {w + a} remaining always between 
Um := min y; and uy := max u;, see e.g. (4.9) or (4.12). Similarly, w takes values in the 
sets {initial values of w} and {u +a}. We conclude that, at every time-step Tk, for almost 
every (x,t) E€ R x [0, 7%) u(x,t) belongs to 


| U U Uu + ka} U fw; + ka} Af lum um] =: LM (u) 


which does not depend on the time-step 7, and necessarily consists of a finite number of 
points. In particular 


OES min dje- y 


rye IM(u) e+ y) >0. (5.8) 


If two discontinuity lines of u cross each other, then the (lower) one must have slope 1/2 
and the (upper) one slope 1 (see Figure 11-right). We denote by uw, u, and u* the values 
of u on the left, on the right and between the two curves, respectively. Such setting tell 
us that the first jump between u, and u* implies a jump for the corresponding pair (u, w) 
between two points belonging to the same boundary of £, and in particular that w changes 
accordingly to u; whereas the second one, from u* to u;, implies a jump for (u, w) between 
two points with the same ordinate w (see Figure 11-left). From this we conclude that u* can 
not be between u and u,. Also note that the Riemann problem solved at the intersecting 
point results in a single discontinuity wave between u; and u, with either slope 1/2 or 1, 
depending on the mutual position of u, and u; and on their distance. In particular, by 
what said above, we have that the total variation of the solution u, in x for t fixed, after 
the intersection point decreases at least by 26, with 6 as in (5.8). Indeed, before it was 
|u, — u*| + |u* — u| and after it is |u, — u|, being |u, — u*| + |u* — u| — Ju, — u| > 26 > 0. 
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Figure 11: Explaining figure for Lemma 5.2. We see how u* cannot be between uw and ur so after the 
interaction the total variation of u decreases. Moreover we generate a 0-slope discontinuity in w as w* # wr. 


Since the total variation of u at time t = 0 is finite and it is non increasing in time (see 
(5.3)) then we cannot have an infinite number of such interactions. 

We finally note that the corresponding jump of w before the intersection was a jump 
between two values w, and w; and after becomes a jump between w, and a new intermediate 
value w* (see Figure 11). In particular new 0-slope discontinuity in w may arise. 


Lemma 5.3. The number of discontinuity lines of u intersecting the discontinuity lines of 
w is finite. 


PROOF. Since the discontinuity lines for w with slope 1/2 are also discontinuity lines for u, 
we restrict ourselves to the discontinuity lines for w with slope 0. 

The discontinuities of w with slope 0 starting at time t = 0 are finite and a new discon- 
tinuity with slope 0 for w may arise only by the interactions of two discontinuity lines for 
u, as in Lemma 5.2. Hence the total number of discontinuities of w with slope 0 that may 
be generated during the whole process is a-priori bounded. Now when a discontinuity line 
of u crosses a 0-slope discontinuity of w, we may generate a new discontinuity of u but the 
number of 0-ones for w does not increase, see the example in Figure 12. 

Since the discontinuity lines for u have positive slope, it is not possible that one of them, 
say rı, intersects a vertical discontinuity line for w that has already generated discontinuity 
lines from which, after a finite number of new intersections, rı itself is generated. 

Concluding, since at time 0 the number of discontinuities for u is finite, then there is a 
finite number of intersection with w-discontinuity of 0 slope. 


Now, we prove the existence of a solution for the Cauchy Problem (3.1), with general 
initial conditions ug, wo. 


Theorem 5.4. Consider (3.1) with T > 0 and uo, wo € BV(R)NL1(R), satisfying |uo(x) — 
wo(x)| < a for almost every x. Then there exists a couple (u, w) € L'(Rr) weak solution in 
the sense of Definition 3.1. 
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Figure 12: Explaining figure for Lemma 5.3. We see how the number of discontinuities of u can increase 
while the total variation Var(u(-,t)) remains constant in time. 


PROOF. Since uo,wo € BV(R) then they are bounded (|uo(x)|, |wo(a)| < M a.e. x), their 
limits at infinity exist and such limits are null as uo, wo are also in L1(R). Moreover, there 


is a sequences {ul}, and fo, of piece-wise constant functions approximating ug and 
wo respectively in L'(R) with 


Var(u\”) <Var(uo) and Var(w\”) <Var(wo) Vn EN. (5.9) 
0 0 


Let (Un,Wn) be the weak solution to the problem with initial data ul” and w? 


Lemma 5.1. Since uy” and wh”), in order to approximate up and wọ in L1(R), must vanish 
at infinity (actually, being the number of pieces finite, their last interval of constancy is a 
semi-line where they are null, they vanish outside a compact set), hence it is easy to see 
that also u,(-,t) and w,(-,t) have the same property, for every fixed t. Applying standard 
procedures for the wave-front tracking method, (see [6], [2]) we get the existence of a function 


u € L! (Rr), continuous from [0,7] into Z7,,.(R), such that (up to a subsequence) 


loc 


) as in 


Un — u, in L'(Rr), as n > oo. (5.10) 


By the same argument, there exists w € L'(Rr), continuous from [0,7] into LŁ.(R), such 


that (up to a subsequence) wn > w in L'(Rr). 

We have to prove that u and w satisfy the weak formulation of the PDE (3.2) and the 
hysteresis relations (3.3), (3.4). 

Since (3.2) is solved by every (un, Wn) then it is also true for the limit because of the Lt 
convergence on Rr. Moreover wy, = F(Un, Wo) for a.e. x, consequently |Wn — Un| < a a.e. 
and hence (3.3) comes from a.e. point-wise convergence. 

It remains (3.4). Because of Proposition 1.3, since w, = F(un), we have that for all n 
and almost every x and t 


eae a ae 
J (a9) 


(0,t) 


(y (2, ») kow (5.11) 
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Here with (2(2,-)) we denote the measure associated to the distribution w/,(x,-), with 
fixed x seen as function of time, and with tilde the right continuous with respect to time 
representative of un and wn. 

Now we fix t € (0,7) such that (5.11) holds for a.e. a. Since for every fixed zx, 
og x,t) ~ 
Wn(Z, “Dio = D a } tn (x) Up at 1,ti) then 


[aena (Free ))= DESS BD (2) — w(e)). 


(0,t) 


Now, if we exclude the 0 measure set of all the the points x given by the union of the 
discontinuity points of wo”) with the points x such that (x,t) lies on a discontinuity line of 
wn then it holds that 


w(x) = lim w,(x,s) = w? (x) 
s—0+F 


and 


wo) (x) = lim w,(x,s) = wlz, t). 
soot 


So by the inequality 


N-1 


Yo ataata) - a) > 


tL 


and by integrating (5.11) over R we deduce 
ae fe fh 9) OWn 
ot 


~~ eh 
2 f (wn(x, t)? — wi (a, t)?) de + a] =" 


(R x (0,t)). (5.12) 


Here with a we denote the measure on R x (0, t) associated to the distributional derivative 
of w, with respect to t seen as function of both x and t. It exists since w, is a function 
of bounded variation of two variables (finite number of constant pieces), moreover it is 
Own — £1 Q (2(z,-)), see [1]. 

Now for every ¢ smooth with compact support in R x (0,t) because of (3.2) we have 


Oun | OWn Oun 7 
(= ðt 2.6) =0 


in a distributional sense and, since Up, Wp are of bounded variation as functions of two 
variables, we can see this relation as an equality between the measures associated to the 
distributions i.e. 

ðw, Oun Oun 


ae BE o 8:18) 
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as measures on R x (0,t). We then have 


Te, i a(t) aa wen a(t) 
=- [7 faena (Frey) æ- ff ien a (Fee) at Oa 
( (0,t) 


25 | * (un(2, t)? = ul (2)?) de — 5 | (tn (+00, t)? — un(=00,t)”) dt 
(0,t) 


Here we have applied similar steps as above (between (5.11)—(5.12)), also exploiting the fact 
that & for fixed t is left continuous with respect to x and that un(+o0, t) = un(—oo, t) = 0. 
Then from (5.12) and (5.14) we conclude that for all n € N and ae. t € (0,T) 


Le 2 p a2 N 2 p a2 JWn 

5 (Un(x,t)“ — ug (£Y) dx + 5 (wnx, t) — wy (£Y) dx +a BE (R x (0,t)) < 0. 
> = (5.15) 

We know by construction and by the wave-front tracking method we have applied (see [6], 


[2]) that for every fixed t the sequences u,,(-,t), Wn(-,t), ul” and wà)” converge respectively 


to u(-,t), w(-,t), uo and wo in L'(R), and since they are all equibounded in L®(R) we get 
the convergence also in L*(R). This in particular means that 


(n) 
0 


0 < lim sup (R x (0,t)) < C < +00, 


n—> o0 


| wn 


o 


and so, up to a subsequence, there exists a measure weak star limit of the sequence Ow, /Ot 
which must then coincide with Ow/dt by the convergence of wp to w in L'(Rr). Finally, 
thanks to the lower-semicontinuity 


Wn 


(R x (0,4)), 


n—> o0 


5 (R x (0,t)) < liminf 


taking the liminf in (5.15), we infer our desired condition (3.4). 


Remark 5.5. We note that in [21] similar arguments, as in the proof above, are used but for 
approximating H+ functions, generated by a time-discretization. Here instead the functions 
are much less regular, the discretization is in space for the initial data, and the derivatives 
are just measures. 
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6. Uniqueness 


Let us recall that the couple (u, w) weak solution is an entropy solution if it satisfies the 
following inequality 


L E Fe E 0: (6.1) 
Rr 


for every ọ non negative, smooth and with compact support in R x (0,T) and for every 
couple (k, È) € £. 

Now we show how in our problem with hysteresis this entropy condition can be character- 
ized in the case of piecewise constant solutions (see e.g. [2] for the case without hysteresis). 


Proposition 6.1. Suppose that the couple of functions u,w: R x [0, +00) > R, satisfying 
ju — w| < a, consists of two constant states (u, w1) # (ur, wr) separated by a line of dis- 
continuity (s(t), t). Then (u, w) is an entropy solution if and only if it is constructed as in 
Section 4 i.e. the following conditions are verified: 


i) s'(t) =: A is constant and satisfies the generalized Rankine — Hugoniot condition (3.5); 

) A is either 0,1 or 1/2; 

iii) if uw =u, thenX=0 i.e. we have a discontinuity only on w with slope 0; 

) if w, = w then X = 1 so the jump in not on the boundary of the hysteresis region and 

w remains constant; 

v) ifu, > u and w, 4 w then A = 1/2 and w, = u, +a, w= u +a, i.e. the jump occurs 
in both u and w and the couple is on the upper boundary of the hysteresis region; 

vi) ifu, < u and w, # w, then A = 1/2 and w, = ur — a, w, = u — a, i.e. the jump occurs 
in both u and w and the couple is on the lower boundary of the hysteresis region; 


PROOF. First notice that integrating by parts (6.1) with (u, w) piece-wise constant we can 
rewrite such entropy condition as follows 


s (E) [lur — kl + [wr — Bl) — (lea — kl + [wr — &I)] — [lu — k| — ju — kl] > 0. (6.2) 


So the couple (u, w) is an entropy solution if and only if (6.2) holds for every (k, Ẹ) € £. 
( == ) Suppose that the pair (u, w) is an entropy solution. Then, for any ur, ur Wr, wi, we 
can always find a pair (k,k) € £ with k, k large enough such that (6.2) reads as follows 


s'(t)[uy + wi — ur — Wy] — [u — ur] > 0 (6.3) 
and, by inserting in (6.2) their opposite values, it becomes 
—s'(t)[u + wi — ur — wy] + [ui — uy] > 0. (6.4) 
We then infer that s’(t) has to satisfy the Rankine-Hugoniot condition (3.5) that is 


[uy — u, + Ww — wrs (t) = wy — Uy, (6.5) 
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hence s/(t) =: À has to be constant, so we proved 7). Notice that u; — u, + w, — w, cannot 
be equal to 0 otherwise from (6.5) we would get also u; — u, = 0 and so w, — w, = 0. This 
means uj = Ur, Wy = wy which we excluded by hypothesis. 

We can now rewrite (6.5) as follows 


Alw, — wy] + (A — 1) (u — u) = 0. (6.6) 


So if u = u, then A = 0, and the left-hand-side of (6.2) becomes 0, independently from 
(k,k). So the couple (u, w) is still an entropy solution and we imply case iii). 

If w, = wr, then A = 1 so we infer case iv). Notice that even in this case (6.2) is trivially 
satisfied. 

Suppose now instead w, Æ w; and up Æ w and recall that u; — ur + w — w, is either 
positive or negative. Keeping this in mind and recalling (6.5) we write (6.2) as follows 


(u — ur)([wr — k| — [wi — Kl) + (wr — wr) (lur — k| — ta = Fl) 
Ul — Up + W — Wp 


> 0. (6.7) 


Now, for (k, k) € R? we study the function 


A 


h(k, k) := (u — Ur) (ler — k| — jwi — bl) + (wy = wr) (leer = k| — [u — kl), 


which is the numerator appearing in (6.7). We can split the plane R? in 9 regions, rep- 
resenting the cases when k is above, below or between wu), up combined with the cases k 
above, below or between wy, wp, see Figure 13. It is easy to prove that h is constant in 
Rı, R3, Ry and Rg as in these parts both k and k cancel out. In particular, in R, it is 
h = 2(wm — u,)(w; — w,) which is strictly positive if the jump between u; and u, has the 
same monotonicity of the jump between w, and w; and strictly negative otherwise. In Ry 
instead h = —2(u; — u,)(w; — wr), and in Rs and Rzy it holds h = 0. It is easy so see then 
that in the other regions h is affine, connecting continuously the constant parts, so we can 
deduce that it has the same constant sign in Rə and R4 and the opposite one in Rg and Rs. 
In R; we have a change of sign on the segment AB connecting (min(w, ur), min(wy, w,)) to 
(max(u;, ur), max(w7, w,)). 

With this in mind we then conclude that if the jump in w has opposite sign with respect 
to the jump in u then this violates our condition (6.2). This is because both the pairs 
(ur, wr) and (w, w) are contained in £ and, see Figure 13, we can find multiple pairs in 
L (eg. (uz, w), (ur, wr) themselves) such that h is either positive or negative. Since the 
denominator u; — u, + w; — w, # 0 has fixed sign this means that the inequality (6.2) cannot 
be true for all (k,k) € £. 

So we can only have jumps in w with the same sign of jumps in u. If u; > u, then the 
denominator u; — Ur + w; — wr is positive, hence it must be h > 0 for all (k, k) € L. But 
the only possibility is that the points A = (up, wr), B = (w,w,) (see Figure 13) both belong 
to the lower boundary of £. Otherwise at least one of them has a neighborhood containing 
points (k, k) in £ where h is negative. So positive jumps occur only on the lower boundary 
of L. Similarly, if u; < u,, then the points must lie on the upper boundary of £. So negative 
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Figure 13: The sign of h, when the jump between u, and ur has the same monotonicity as the jump between 
w and wp. Here A = (min(u, ur), min(wr, w1)) and B = (max(uz, ur), max(w,;, wi)). Note that if the jump 
between u; and u, has different monotonicity to the jump between w, and wr, then the picture remains the 
same but with opposite signs, in particular, the segment AB is the same. In this case the pairs (uw, wr), 
(ur, wr) lie on the highlighted points. 


jumps occur only on the upper boundary of £. Also note that, in both cases, u and w must 
have a jump of the same amplitude. 

So we checked that for w 4 ur and w, 4 wr, (6.2) implies A = 1/2 (see (6.5)) and either 
case v) or vi). 

As there are no other possible choices of (uw, w7), (up, wp) we also infer ii). 
( <= ) Now suppose 7),..., vi) to hold. If u, = u, then we are in case iii) with A = 0. As 
already said, the left-hand-side of (6.2) becomes 0, hence the inequality is trivially satisfies 
for every (k,k) € £. 

If w; = w, then we are in case iv) with \ = 1 and again (6.2) is trivially satisfies. 

Suppose instead A = 1/2, and notice again, because of i), that u; — u, + w, — w, # 0 
otherwise we would get u; = u, and w, = w, which we excluded by hypothesis. This allows 
us to rewrite (6.2) as (6.7). Applying the same reasoning as for the previous implication, it 
can be verified that for both the cases v) and vi), (6.7) holds true for every (k, Ẹ) € £. 


Thanks to Proposition 6.1, it is easy to see that the sequence of piece-wise constant solutions 
(Un; Wn) of the Cauchy problem with ul and wi” as initial data constructed via the Wave- 
Front Tracking, is entropy admissible, so it satisfies (6.1). Finally, since we have L1(Rr) 
convergence, we can conclude that the solution (u, w) as constructed in the proof of Theorem 
5.4 also satisfies (6.1). As a consequence of the next theorem, we then deduce that such a 
solution is also the only entropy solution. 


Theorem 6.2. Consider the Cauchy problems with initial conditions uj, w} and ug, we re- 
spectively, where uj, wi, € L'(R)OBV (R). Let us denote by (u1, w1) and (u2, w2) two entropy 
solutions to the respective problems. Then it holds that 


+00 


+00 
f (lui = val(æ,t) + lun — walle) da < f (jii — a2] + Jw aede (6.8) 


for almost every t € [0, T]. 
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Proor. Adapting the standard method by Kruzkov [14], in [21] an analogous result is 
proven for the case of delayed-relay hysteresis. The adaptation to our case with Play hys- 
teresis is not difficult. 
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